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Nuclear Forensics and Radiochemistry: Reaction Networks 

Robert S. Rundberg 


Abstract: 

In the intense neutron flux of a nuclear explosion the production of isotopes may 
occur through successive neutron induced reactions. The pathway to these isotopes 
illustrates both the complexity of the problem and the need for high quality nuclear 
data. The growth and decay of radioactive isotopes can follow a similarly complex 
network. The Bateman equation will be described and modified to apply to the 
transmutation of isotopes in a high flux reactor. A alternative model of growth and 
decay, the GD code, that can be applied to fission products will also be described. 
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Reaction Networks 


A nuclear explosion produces an enormous flux 
of neutrons. 

- A neutron yield on the order of a mole (6.02E23). 

- The chain reaction is finished in a short time, of the 
order of tens of nano-seconds. 

- The volume of burning fuel is relatively small. 

Multiple successive neutron induced reactions 
can occur. 

Activation of materials by successive reactions 
are prominent near the fuel. 



Neutron Reactions with Actinides 
(Uranium) 








































Neutron Reactions with Actinides 
(Uranium) 
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Neutron Reactions with Actinides 
(Plutonium) 




































Neutron Reactions with the Actinides 
(Americium) 










































Actinide Isotopes in Nuclear Debris 


• The equations for the successive reactions are 
solved numerically in a code using the Runga- 
Kutta method. 

• The isotopes produced give the neutron 
fluence and spectrum that the original 
material was exposed to. 



The Advantage of Isotope Ratios 

If we assume no burnup the first order equation is 


dNA—i 
dt 


= <71 N A (j). 


Assume constant flux, the solution is 


N A -i = (J\N A (j)t . 


The second order reaction is then 

dN ^ t ~ 2 = ^2 [<JlN A (j)t] (/). 


(1) 

(2) 

( 3 ) 


The solution is 


N A -2 = CTl(T2N A [ (j) 2 tdt 
Jo 

= cr^NA^t ) 2 
2 


( 4 ) 


Na-2 
N a _i 



Therefore the atom ratio is 


(5) 



Growth and Decay 


The growth and decay of radioactive isotopes 
are similar to the transmutation of isotopes in 
high flux environment. 

The Bateman equations are useful. But 
restricted to specific initial conditions. 

The Bateman equations can be modified to 
include transmutation. 
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23S U-series 


23S U-series 


232 Th-series 
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Bateman Equations 


Start with a simple parent/daughter growth and decay. The equation for the 
daughter is, 

^ = \\Ni - X 2 N 2 
at 


or 


dN 2 

dt 


+ X 2 N 2 - XiN°e~ Xl = 0. 


( 1 ) 


The solution of this linear differential equation of the first order may be obtained 
by standard methods and gives 


Ai 


A2 — A] 


9 _j_ 


-AT? (, 


N 2 = 


( 2 ) 




Bateman Equations 


Consider the grandaughter. The equation for its growth and decay is, 


~ = a 2 at 2 - x 3 n 3 . 


(3) 


Eq. 3 is analgous to Eq. 1, but the solution calls for more labor, because AT 2 
is a much more complicated function than N\. The great grandaughter is even 
more complicated. Fortunately, H. Bateman has given the solution for a chain 
of n members with the special assumption that at t = 0 the parent substance 
alone is present, that is, = ■ • • = N® = 0. This solution is 


N n = C ie ~ Xlt + C 2 e~ X2t +■■■ C n e~ Xnt , 


(A 2 - Al)(A 3 - Al) • * • - Al) 




JV^and so on. 


( 4 ) 




Bateman Applied to Transmutation 


Successive neutron reactions in a high flux, such as, a high flux reactor can be 
solved using the Bateman equations, as well. The rate of disappearence of an 
isotope in a neutron flux is 


—— = (A + nva)N = AN. (5) 

at 

Consider a parent daughter pair. The parent disappears by both transmutation 
and decay. But the daughter grows by decay of the parent only and disappears 
by both processes. In general notation, 

= XiK - A i+1 N i+1 . 

We replace A* by a modified decay constant A* = A* + nz/cr*, were only the 
decay constant of transmuation term that lead to next progeny in the chain is 
used. 



Bateman Applied to Transmutation 


With this nomenclature the Bateman equation becomes, 
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( 6 ) 




An Example 


As an illustration, we compute the amount of 3.15-d 199 Au formed by two 
successive (n, y) reactions when 1 g 197 Au is exposed for 30 h in a neutron 
flux of 1 x 10 14 cm' 2 s" 1 . The chain of reactions is 


7 Au- 


8 Au - 


<1/2-2.70 d 

We use (5-12) for this three-membered chain: 
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The numerical values to be substituted are 

t = 1.08 x 10 5 s, 
nv = 10 14 cm” 2 s“', 
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A 1*97 
A 198 

A 198 


= 9.9 x 10 -23 cm 2 , 
= 2.5 x 10 -20 cm 2 . 


6.02 x 10 23 
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and 



Using these values, we get 

( _ -0.00107 

5.5 x 10~* x 2.55 x I0~* 

e -0.594 e -0.275 ^ 

+ 5.5x 10- 4 x 2.95 x 10 -6- 2.55 x 10' 6 x 2.95 x 1(TV 

= 7.55 x 10 7 (7.J2 x 10 ,o + 3.40 x 10'°- 1.01 x 10") = 3.2 x 10 17 . 

The disintegration rate of ,99 Au at the end of the irradiation is A 199^,99 = 
0.82 x 10 12 s~ ! . For comparison we compute the disintegration rate of ,98 Au 
in the sample [again from (5-12) for a two-membered chain]: 

/ g -A 197» \ 
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= 9.06 x 10 7 
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Thus about 10 percent of the radioactive disintegrations in the sample occur in 
199 Au. 



Fission Does Not Always Meet the 
Initial Conditions 






















The GD Code 


The general differential equations for radioactive decay and growth 
dNi 

-jj-=a n Ni +n(t) 

— tt = a>2iNi -f Q 22 N 2 +r2(t) (1) 


— a n\N\ +-h a nnNn +^*n(^) 

where Ni(t) is the number of atoms of nuclide i existing at time £, a,ij , i > j are 
constants, and n J s are the rates of formation from sources other than by decay 
of isotopes i. A specialized form t\(£) = yif(t) is used in the Los Alamos code 
GD, so that the rates have time independent ratios to each other. The y^s are 
called fractional independent yields. 


GD Code Continued 


where A is the matrix of a^’s, N and Y are column vectors and f(t) is a 
scalar. Equation (2) has the solution 

Is (3) 


N(t) = e<‘- T ) A N(r) + J e^ A f{t)Yds 
There is a non singular matrix P such that 

P -1 AP = J 

is diagonal. Th GD code finds P. The solution is then 

N(i) = P diag [e ail (*- T ),..., e a ""( t - T >] P _1 N(t) 


(4) 
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The GD Code Solution 


Define Fi(t) by 

Fi(t — t) = j l e a ‘ dt ~ s) f(s)ds (6) 

then the solution is 


N (t) = 

P diag , e 0 ""( t -' r >] P _1 N(t) 

+ 

P diag [Pi(t - t), .. .,F„(t - t)] P X Y (7) 


The function f(s) is taken to be piecewise constant with fj on time intervals 
Uj to Wj . With this specialization it follows that 


m -r) =—Y,fj ( e ~ atiWl - e ~ auUi ) • 

(in 


( 8 ) 



